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Abstract. In this survey a novel system, called the oscillator system, con- 
sisting of order of p"^ functions (signals) on the finite field Fp, is described and 
studied. The new functions are proved to satisfy good auto-correlation, cross- 
correlation and low peak-to-average power ratio properties. Moreover, the 
oscillator system is closed under the operation of discrete Fourier transform. 
Applications of the oscillator system for discrete radar and digital communi- 
cation theory are explained. Finally, an explicit algorithm to construct the 
oscillator system is presented. 



1. Introduction 

One-dimensional analog signals are complex valued functions on the real line M. 
In the same spirit, one-dimensional digital signals, also called sequences, might be 
considered as complex valued functions on the finite line Fp, i.e., the finite field 
with p elements, where p is an odd prime. In both situations the parameter of 
the line is denoted by t and is referred to as time. In this survey, we will consider 
digital signals only, which will be simply referred to as signals. The space of signals 
TL =C(Fp) is a Hilbert space with the Hermitian product given by 



t£¥p 

A central problem is to construct interesting and useful systems of signals. Given 
a system S, there are various desired properties which appear in the engineering 
wish list. For example, in various situations [ll[2] one requires that the signals will 
be weakly correlated, i.e., that for every (j) ip £ & 

1(0,^)1 «1- 

This property is trivially satisfied if 6 is an orthonormal basis. Such a system 
cannot consist of more than dimTi signals, however, for certain applications, e.g., 
CDMA (Code Division Multiple Access) [3] a larger number of signals is desired, 
in that case the orthogonality condition is relaxed. 

During the transmission process, a signal ip might be distorted in various ways. 
Two basic types of distortions are time shift ip(t) i—s- L^(^(<) — (p{t + r) and phase 
shift ip{t) I—)- Mwip{t) = e p ^ where t,w & Fp. The first type appears in 

asynchronous communication and the second type is a Doppler effect due to relative 
velocity between the transmitting and receiving antennas. In conclusion, a general 
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distortion is of the type ip i-+ MmLr^p, suggesting that for every ^ e S it is 
natural to require [2] the following stronger condition 

\{(j)Mw\-r^p)\ < 1- 

Due to technical restrictions in the transmission process, signals are sometimes 
required to admit low peak-to-average power ratio [4], i.e., that for every Lp & & 
with \\p\\2 = 1 

WLKK{\(f{t)\ : t e Fp} < 1. 

Finally, several schemes for digital communication require that the above properties 
will continue to hold also if we replace signals from 6 by their Fourier transform. 

In this survey we demonstrate a construction of a novel system of (unit) signals 
So, consisting of order oi signals, called the oscillator system. These signals 
constitute, in an appropriate formal sense, a finite analogue for the eigenfunctions 
of the harmonic oscillator in the real setting and, in accordance, they share many 
of the nice properties of the latter class. In particular, the system &o satisfies the 
following properties 

(1) Auto- correlation (ambiguity Junction). For every ip G So we have 

, , r 1 if (r,-u;) 0, , , 

i(^,M.u^>i = | if (;,;,)Vo. (1-1) 

(2) Cross-correlation (cross-ambiguity function). For every (p ^ ip ^ So we 
have 

|(<^,M^U^)| < A, (1.2) 

for every r, w e Fp. 

(3) Supremum. For every signal ip S So we have 

2 

max{|v3(t)| -.t e¥„} < —. 

Vp 

(4) Fourier invariance. For every signal p € So its Fourier transform p is (up 
to multiplication by a unitary scalar) also in So- 

The oscillator system can be extended to a much larger system S^, consisting 
of order of p^ signals if one is willing to compromise Properties 1 and 2 for a 
weaker condition. The extended system consists of all signals of the form M^L^ip 
for T, w g Fp and p G So- It is not hard to show that # (S^) = • # (So) ~ p^- 
As a consequence of (|l.ip and (|1.2p for every 7^ e S^ we have 

l(.,0)l<^- 

The characterization and construction of the oscillator system is representation 
theoretic and we devote the rest of the survey to an intuitive explanation of the 
main underlying ideas- As a suggestive model example we explain first the con- 
struction of the well known system of chirp (Heisenberg) signals, deliberately taking 
a representation theoretic point of view (see [21 [5] for a more comprehensive treat- 
ment). 
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2. Model example (Heisenberg system) 

Let us denote hy tp : ¥p the character ■ip{t) = e~*. We consider 

the pair of orthonormal bases A = {6a ■ a G Fp} and = {■^^ : a G Fp}, where 

2.1. Characterization of the bases A and A^. Let L : 7i — !■ 7i be the time shift 

operator L(p{t) = (p{t + l). This operator is unitary and it induces a homomorphism 
of groups L : Fp ^ U{H) given by LTfit) = ^{t + r) for any r € Fp. 

Elements of the basis A^ are character vectors with respect to the action L, 
i.e., \-Ttpa ~ V'(^''')V'a ^ ^p- same fashion, the basis A consists of 

character vectors with respect to the homomorphism M : Fp ^ [/ (TL) generated by 
the phase shift operator M^p{t) = il){t)(p{t). 

2.2. The Heisenberg representation. The homomorphisms L and M can be 

combined into a single map n : ¥p x ¥p ^ U {Ti.) which sends a pair (t, w) to the 
unitary operator tt{t,lo) = ip (^— ^tw) M^, o L^. The plane Fp x Fp is called the time- 
frequency plane and will be denoted by V. The map tt is not an homomorphism 
since, in general, the operators and do not commute. This deficiency can 
be corrected if we consider the group H = V x ¥p with multiplication given by 

(r, w, z) ■ (r', w' , z') = (t + r', w + w', z + z' + \{tw' — t'w)). 

The map tt extends to a homomorphism tt : H U{7{) given by 

7r(T, w,z) = -ip {—-^TW + z) Mm o Lt- 

The group H is called the Heisenberg group and the homomorphism n is called the 
Heisenberg representation. 

2.3. Maximal commutative subgroups. The Heisenberg group is no longer 
commutative, however, it contains various commutative subgroups which can be 
easily described. To every line L C V, which pass through the origin, one can 
associate a maximal commutative subgroup Al = {(/, 0) € V x ¥p : I G L}. It will 
be convenient to identify the subgroup Al with the line L. 

2.4. Bases associated with lines. Restricting the Heisenberg representation tt 
to a subgroup L yields a decomposition of the Hilbert space H into a direct sum 
of one-dimensional subspaces H =0^^, where x runs in the set of (complex 

X 

valued) characters of the group L. The subspace consists of vectors gH such 
that 7r(Z)(/3 = xO-)^- In other words, the space Ti.^ consists of common eigenvectors 
with respect to the commutative system of unitary operators {t^{1)}i^^ such that 
the operator tt (l) has eigenvalue x (0- 

Choosing a unit vector (p^ G Tiy- for every x G wc obtain an orthonormal 
basis Bl = {f^ : x G -^^}- In particular, A^ and A arc recovered as the bases 
associated with the lines T = {(r, 0) : r e Fp} and W = {(0, w) : w G Fp} respec- 
tively. For a general L the signals in Bl arc certain kind of chirps. Concluding, 
we associated with every line L C V a.n orthonormal basis Bl, and overall we 
constructed a system of signals consisting of a union of orthonormal bases 

eH={v^BL:LcV}. 

For obvious reasons, the system &h will be called the Heisenberg system. 
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2.5. Properties of the Heisenberg system. It will be convenient to introduce 
the following general notion. Given two signals G 7i, their matrix coefRcient is 
the function ra^^^p : H ^ C given by m^^^{h) — ((f),iT{h)(p). In coordinates, if we 
write h — (r, w, z) then m^ yj{h) = ip {—^tw + z) (0, M^, o L^iys). When (f> — (p the 
function m^p ,^ is called the ambiguity function of the vector ip and is denoted by 

The system &h consists of p + 1 orthonormal base^, altogether p{p + 1) signals 
and it satisfies the following properties [SJ [5] 

(1) Auto-correlation. For every signal ip e Bl the function is the charac- 
teristic fmiction of the line L, i.e., 



iv)\ 



0, V i L, 

1, w e i. 



(2) Cross- correlation. For every £ Bl and ip G Bm where L ^ M we have 

for every v £ V. If L = M thenjTOi^^^l is the characteristic function of some 
translation of the line L. 

(3) Supremum. A signal ip S &h is a unimodular function, i.e., \'p3{t)\ — 
for every t £ ¥p, in particular we have 



max{|(^(i)| : < e Fp} = — < 1. 



1 

VP 



Remark 1. Note the main differences between the Heisenberg and the oscillator 
systems. The oscillator system consists of order ofp^ signals, while the Heisenberg 
system consists of order of p^ signals. Signals in the oscillator system admit an 
ambiguity Junction concentrated at Q £ V (thumbtack pattern) while signals in the 
Heisenberg system admit ambiguity function concentrated on a line. 



3. The oscillator system 

Reflecting back on the Heisenberg system we see that each vector ip £ &h 
is characterized in terms of action of the additive group Ga — Fj,. Roughly, in 
comparison, each vector in the oscillator system is characterized in terms of action 
of the multiplicative group Gm = ■ Our next goal is to explain the last assertion. 
We begin by giving a model example. 

Given a multiplicative characteiu x ■ Gm — > C ^ , we define a vector x e 7i by 



xit) 



0, t = 0. 



We consider the system Bstd = {x • X S G^, X 7^ l}> where G^ is the dual group 
of characters. 



^Note that p + 1 is the number of Unes in V. 

multipUcative character is a function x '■ Gm —* which satisfies xi^v) = x(^)x{y) foi' 
every x,y e Gm- 
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3.1. Characterizing the system Bstd- For each element a G Gm let : H ^ H 

be the unitary operator acting by scahng p^ip{t) — ip{at). This collection of opera- 
tors form a homomorphism p : Gm C/(7i). 

Elements of Bstd are character vectors with respect to p, i.e., the vector x satisfies 
Pa (x) — x{'^)x for every a £ Gm- In more conceptual terms, the action p yields 
a decomposition of the Hilbert space H into character spaces 7i = ® TCx , where x 
runs in the group G^. The system Bstd consists of a representative unit vector for 
each space H^, x 7^ 1- 

3.2. The Weil representation. We would like to generalize the system Bstd in 
a similar fashion to the way we generalized the bases A and in the Heisenberg 
setting. In order to this we need to introduce several auxiliary operators. 

Let p^ : H H, a £ ¥p , he the operators acting by Pa^{t) = cr{a)ip{a~-^t) 
(scaling) , where a is the unique quadratic character of ¥p , let pj^ : H ^ H be 
the operator acting by prpip{t) = ip{t'^)(p(t) (quadratic modulation), and finally let 
Pg : Ti. —> Ti. he the operator of Fourier transform 

VP seFp 

where i/ is a normalization constant [5]. The operators PaT Pt ^^'^ Ps unitary. 
Let us consider the subgroup of unitary operators generated by Pa, Ps a-^i^ Pt- This 
group turns out to be isomorphic to the finite group Sp = 5*^2 (Fp), therefore we 
obtained a homomorphism p : Sp t/(7i). The representation p is called the Weil 
representation [7| and it will play a prominent role in this survey. 

3.3. Systems associated with maximal (split) tori. The group Sp consists of 
various types of commutative subgroups. We will be interested in maximal diago- 
nalizable commutative subgroups. A subgroup of this type is called maximal split 
torus. The standard example is the subgroup consisting of all diagonal matrices 

which is called the standard torus. The restriction of the Weil representation to a 
split torus T C Sp yields a decomposition of the Hilbert space TL into a direct sum 
of character spaces Ti Hx , where x runs in the set of characters . Choosing 
a unit vector ip^ G Ti.^ for every x we obtain a collection of orthonormal vectors 
St = {<y3^ : X G T^, X 7^ Overall, we constructed a system 

e'o= {(p e Bt ■■ T C Sp split} , 

which will be referred to as the split oscillator system. We note that our initial 
system Bstd is recovered as Bstd — Ba- 

3.4. Systems associated with meuximal (non-split) tori. From the point of 
view of this survey, the most interesting maximal commutative subgroups in Sp 
are those which are diagonalizable over an extension field rather than over the base 
field Fp. A subgroup of this type is called maximal non-split torus. It might be 
suggestive to first explain the analogue notion in the more familiar setting of the 
field M. Here, the standard example of a maximal non-split torus is the circle group 
50(2) C 5L2(R)- Indeed, it is a maximal commutative subgroup which becomes 
diagonalizable when considered over the extension field C of complex numbers. 
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The above analogy suggests a way to construct examples of maximal non-split 
tori in the finite field setting as well. Let us assume for simplicity that —1 does not 
admit a square root in Fp. The group Sp acts naturally on the plane V — ¥p x ¥p. 
Consider the symmetric bilinear form BonV given by 



An example of maximal non-split torus is the subgroup Tns C Sp consisting of 
all elements g £ Sp preserving the form B, i.e., g € Tns if and only if B{gu,gv) = 
B{u, v) for every u,v & V . 

In the same fashion like in the split case, restricting the Weil representation 
to a non-split torus T yields a decomposition into character spaces Ti —^Ti-^- 
Choosing a unit vector ip-^ g Ti.^ for every x G we obtain an orthonormal basis 
Bt- Overall, we constructed a system of signals 



The system Gq' will be referred to as the non-split oscillator system. The con- 
struction of the system &o = &o U together with the formulation of some of 
its properties are the main contribution of this survey. 

3.5. Behavior under Fourier transform. The oscillator system is closed under 
the operation of Fourier transform, i.e., for every if G &o we have if € &o- The 
Fourier transform on the space C (Fp) appears as a specific operator p (w) in the 
Weil representation, where 



Given a signal ip G Bt C Sq, its Fourier transform = p (w) ip is, up to a unitary 
scalar, a signal in Bt' where T' — wTw^^ . In fact, &o is closed under all the 
operators in the Weil representation! Given an element g € Sp and a signal tp € Bt 
we have, up to a unitary scalar, that p (g) p (£ Bt' , where T' — gTg^^. 

In addition, the Weyl element w is an element in some maximal torus (the 
split type of depends on the characteristic p of the field) and as a result signals 
(p e Bt„ are, in particular, eigenvectors of the Fourier transform. As a consequences 
a signal (p £ Bt„ and its Fourier transform ip differ by a unitary constant, therefore 
are practically the "same" for all essential matters. 

These properties might be relevant for applications to OFDM (Orthogonal Fre- 
quency Division Multiplexing) 8 where one requires good properties both from the 
signal and its Fourier transform. 

3.6. Relation to the harmonic oscillator. Here we give the explanation why 
functions in the non-split oscillator system 6^" constitute a finite analogue of the 
eigenfunctions of the harmonic oscillator in the real setting. The Weil representation 
establishes the dictionary between these two, seemingly, unrelated objects. The 
argument works as follows. 

The one-dimensional harmonic oscillator is given by the differential operator 
D = d'^-t^. The operator D can be exponentiated to give a unitary representation 
of the circle group p : SO (2,R) — > U (^^(M)) where p{0) = e^^^. Eigenfunctions 
of D are naturally identified with character vectors with respect to p. The crucial 
point is that p is the restriction of the Weil representation of SL2 (M) to the maximal 
non-split torus SO (2,R) C SL2 (M). 



B{{t,w),(t',w')) ^tt' + ww'. 



{peBt-.T C Sp non-spfit} . 
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Summarizing, the cigcnfunctions of the harmonic oscillator and functions in ©q* 
are governed by the same mechanism, namely both are character vectors with re- 
spect to the restriction of the Weil representation to a maximal non-split torus 
in SL2. The only difference appears to be the field of definition, which for the 
harmonic oscillator is the reals and for the oscillator functions is the finite field. 

4. Applications 

Two applications of the oscillator system will be described. The first application 
is to the theory of discrete radar. The second appfication is to CDMA systems. We 
will give a brief explanation of these problems, while emphasizing the relation to 
the Heisenberg representation. 

4.1. Discrete Radar. The theory of discrete radar is closely related [2] to the 
finite Heisenberg group H. A radar sends a signal Lp{t) and obtains an echo e{t). 
The goal ^ is to reconstruct, in maximal accuracy, the target range and velocity. 
The signal Lp{t) and the echo e{t) are, principally, related by the transformation 

e{t) = e2™V(i + 7-) = M^U(^(t), 

where the time shift t encodes the distance of the target from the radar and the 
phase shift encodes the velocity of the target. Equivalently saying, the transmitted 
signal and the received echo e are related by an action of an element hg G H , 
i.e., e = TT{ho)(p. The problem of discrete radar can be described as follows. Given 
a signal ip and an echo e — Tr{h())ip extract the value of hg. 

It is easy to show that |mi^_e (^)| = l^tp {h ■ ho)\ and it obtains its maximum 
at Hq^. This suggests that a desired signal (p for discrete radar should admit an 
ambiguity function A^p which is highly concentrated around £ H, which is a 
property satisfied by signals in the oscillator system (Property 2). 

Remark 2. It should be noted that the system Go is "large" consisting of aproximately 
signals. This property becomes important in a jamming scenario. 

4.2. Code Division Multiple Access (CDMA). We are considering the follow- 
ing setting. 

• There exists a collection of users i G I, each holding a bit of information 
fei € C (usually bi is taken to be an iV'th root of unity). 

• Each user transmits his bit of information, say, to a central antenna. In 
order to do that, he multiplies his bit bi by a private signal ip^ € Ti. and 
forms a message Ui — biip^. 

• The transmission is carried through a single channel (for example in the 
case of cellular communication the channel is the atmosphere), therefore 
the message received by the antenna is the sum 

i 

The main problem [3] is to extract the individual bits bi from the message u. 
The bit bi can be estimated by calculating the inner product 

The last expression above should be considered as a sum of the information bit 
bi and an additional noise caused by the interference of the other messages. This 
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is the standard scenario also called the Synchronous scenario. In practice, more 
complicated scenarios appear, e.g., asynchronous scenario - in which each message 
Ui is allowed to acquire an arbitrary time shift Ui{t) i— > Ui{t + Ti)^ phase shift 
scenario - in which each message Ui is allowed to acquire an arbitrary phase shift 
Ui{t) ^ e~^^*Ui{t) and probably also a combination of the two where each message 
Ui is allowed to acquire an arbitrary distortion of the form Ui(t) t-^ e~'^^*Ui{t + Ti). 

The previous discussion suggests that what we are seeking for is a large system 
S of signals which will enable a reliable extraction of each bit bi for as many users 
transmitting through the channel simultaneously. 

Definition 1 (Stability conditions). Two unit signals (j) ^ ip are called stably 
cross- correlated if \m^^^{v)\ <ti 1 for every v d V . A unit signal ip is called 
stably auto- correlated if\A^ {v)\ <^ 1, forv ^ 0. A system 6 of signals is called 
a stable system if every signal if €z & is stably auto- correlated and any two different 
signals 4>,(p €i & are stably cross-correlated. 

Formally what we require for CDMA is a stable system 6. Let us explain why 
this corresponds to a reasonable solution to our problem. At a certain time t the 
antenna receives a message 

which is transmitted from a subset of users J C I. Each message Ui, i G J, is of 
the form Ui = bie~"''^^tp^{t + t^) = biTT{hi){p^, where hi G H. In order to extract 
the bit bi we compute the matrix coefficient 

= b^RhiA^. + #(J - {i})o(l), 

where Rh^ is the operator of right translation Rh^Aip.{h) = Aip.{hhi). 

If the cardinality of the set J is not too big then by evaluating r7i<p.^„ aX h — h^^ 
we can reconstruct the bit bi. It follows from (jl.ip and (|1.2p that the oscillator 
system &o can support order of users, enabling reliable reconstruction when 
order of y/p users are transmitting simultaneously. 



Appendix A. Algorithmic construction of the oscillator system 

A.l. Algorithm. We describe an explicit algorithm that generates the oscillator 
system &q associated with the collection of split tori in Sp. 

A. 1.1. Tori. Consider the standard diagonal torus 

A^{{; ;.);<.e.;}. 

Every split torus in Sp is conjugated to the torus A, which means that the 
collection T of all split tori in Sp can be written as 

r = {gAg-^; g e Sp}. 
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A. 1.2. Parametrization. A direct calculation reveals that every torus in T can be 
written as gAg^^ for an element g of the form 

If & = 0, this presentation is unique: In the case 6 7^ 0, an element g represents 
the same torus as g if and only if it is of the form 

'1 & \ / 

c l + bc) [b-'^ 

Let us choose a set of elements of the form (|A.1I) representing each torus in T 
exactly once and denote this set of representative elements by R. . 

A. 1.3. Generators. The group A is a cyclic group and we can find a generator 
for A. This task is simple from the computational perspective, since the group A 
is finite, consisting of p — 1 elements. 

Now, we make the following two observations. First observation is that the 
oscillator basis Ba is the basis of eigenfunctions of the operator p {gA)- 

The second observation is that, other bases in the oscillator system &q can be 
obtained from Ba by applying elements from the set R. More specifically, for a 
torus T of the form T = gAg^^, g E R, we have 

Concluding, we described the (split) oscillator system 

&h={pi9)v-9^R,V^BA}. 

A. 2. Formulas. We are left to explain how to write explicit formulas (matrices) 
for the operators p{g), g G R. 

First, we recall that the group Sp admits a Bruhat decomposition Sp = BUBwB, 
where B is the Borel subgroup consisting of upper triangular matrices in Sp and w 
denotes the Weyl element 

• 1^ 

Furthermore, the Borel subgroup B can be written as a product B = AU ~ UA, 
where A is the standard diagonal torus and U is the standard unipotent group 



Therefore, we can write the Bruhat decomposition also as Sp = UAU U AwU. 

Second, we give an explicit description (which can be easily verified) of operators 
in the Weil representation which are associated with different types of elements in 
Sp. The operators are specified up to a unitary scalar, which is enough for our 
needs. 

• The standard torus A acts by (normalized) scaling: An element 



^0 

acts by 

Sa [/] (t) = a{a)f (a-H) , 
where ct : ^ {±1} is the Legendre character, a{a) = a^(modp). 
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• The subgroup of strictly lower diagonal elements U (Z Sp acts by quadratic 
exponents (chirps): An element 

,u 1, 



u = 



Mulf] (i)=V(-f 



w = 



acts by 

Mu [f] (i) = V-l- 2 ' 

where ^ : Fp — > is the character tjj{t) = e"?"*. 
The Weyl element 

1 
^-1 

acts by discrete Fourier transform 

Using the Bruhat decomposition we conclude that every operator p{g), 
g G Sp, can be written either in the form p [g) = M„ o Sa or in the form 
p {g) = M„2 o SaO F o , where M„, 5a and F are the explicit operators 
above. 

Example 1. For g G R, with 6^0, the Bruhat decomposition of g is given 
explicitly by 

_ ( I 0\ fh \ f 1\ f I 0^ 

and consequently 

p [g) = Mi+bc o Sbo Fo Mb-i . 

b 

For g G R, with c = 0, we have 

'1 0^ 
u 1 

and 

P (g) = Mu. 

A. 3. Pseudocode. Below, is given a pseudo-code description of the construction 
of the (split) oscillator system &q. 

(1) Choose a prime p. 

(2) Compute generator gA for the standard torus A. 

(3) Diagonalize p{gA) and obtain the basis of eigenfunctions Ba- 

(4) For every g G R: 

(5) Compute the operator p (g) as follows: 

(a) Calculate the Bruhat decomposition of g, namely, write g in the form 

g = U2 ■ a ■ w ■ ui or g = u ■ a. 

(b) Calculate the operator p (g), namely, take p (<?) = M„2 o Sa o F o Mu^ 
ov p{g) = M„ o Sa- 

(6) Compute the vectors p{g)<fi, for every (fi G Ba and obtain the system 

BgAg-^- 



Remark 3 (Running time). It is easy to verify that the time complexity of the 
algorithm presented above is 0{p^logp). This is, in fact, an optimal time 
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Remark about field extensions. All the results in this survey were stated for 
the basic finite field ¥p for the reason of making the terminology more accessible. 
However, they are valid for any field extension of the form ¥q with q — p". Complete 
proofs appear in [6 . 
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